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Abstract. The production of charmed particles by E~ of 340 GeV / c momentum was studied in the hyperon 
beam experiment WA89 at the CERN-SPS, using the ^-spectrometer. In two data-taking periods in 1993 
and 1994 an integrated luminosity of 1600 (ib^ 1 on copper and carbon targets was recorded. From the 
reconstruction of 930 ± 90 charm particle decays in 10 decay channels production cross sections for D, D, 
DjTand A+ were determined in the region xf > 0. Assuming an A 1 dependence of the cross section on 
the nucleon number, we calculate a total cc production cross section of o c -(xf > 0) = 5.3 ± 0.4 (stat) ± 
1.0 (syst) +1.0 (H c )/ib per nucleon. The last term is an upper limit on the unknown contribution from 
charmed-strange baryon production. 



2 



1 Introduction 

Nearly 25 years after the discovery of the charm quark, 
charm hadron physics is still a major field of research. 
This results from the difficult identification of charm par- 
ticle decays due to low production cross sections, short 
lifetimes and branching ratios of only a few percent in 
the principal decay modes. Not only the properties of the 
charmed hadrons themselves are the object of intense ex- 
perimental studies. Charm quark production is also of in- 
terest, since it can be used for tests of QCD calculations. 
As the lightest of the "heavy" quarks, the charm quark is 
of special interest as its production can be described using 
perturbative methods and it is experimentally accessible 
with reasonable statistics. 

Usually charm particle production is modeled as a two- 
step process involving different energy scales: one hard 
part, which describes the production of the charm quark 
pair itself, and the soft process of subsequent hadroniza- 
tion. 

The elementary charm quark production process, dom- 
inated at typical fixed target energies by gluon-gluon fu- 
sion, has been calculated in the framework of perturbative 
QCD in leading and next to leading order, the NLO con- 
tributions exceeding the leading ones by a factor of three 
(see, for instance, Q). 

Because of the relatively small charm quark mass, NLO 
calculations still contain considerable uncertainties com- 
ing from uncertainties in the renormalization and factor- 
ization scales, unprecise parton density functions as well 
as the unknown charm quark mass. At present these un- 
certainties are much larger than the experimental errors. 
Therefore measurements of charm production in differ- 
ent beams and at different energies provide a data set on 
which many further developments of QCD-based calcula- 
tions can be tested rigorously. 

In fixed target experiments charm particle production 
has so far been studied mainly with pion and proton beams, 
with beam momenta ranging from 200 to 600 and 200 to 
800 GeV/c , respectively. Only few data on charm pro- 
duction by a kaon beam exist so far. All measured cross 
sections lie within the wide margins of the NLO calcula- 
tions. Many of these experiments measured D/D produc- 
tion cross sections only, and the total charm production 
cross sections have to be extrapolated from these measure- 
ments. 

In this paper we report on the first measurement of 
the total cc cross section in E~ - nucleus interactions at 
y/s = 25GeV/c 2 , based on observed samples of D ± , D°, 
D°, D~, and A+. Results on the charge asymmetries in 



the production of D + , D° , D+ A c and their antiparticles 
have already been published H . 



2 Experimental setup 

Experiment WA89 at the CERN-SPS used a secondary 
beam of created in 450 GeV/c p-Be interactions. The 
detector was designed as a typical forward spectrometer, 
set up in the CERN West area around the /2-magnet, 
and consisted of a target region, a decay region, the fi- 
spectrometer, a ring imaging Cherenkov counter and an 
electromagnetic and a hadronic calorimeter. Figure [l] shows 
the setup in 1993 - the setup in 1994 differed only margi- 
nally. 

The hyperon beam was produced from 450 GeV/c pro- 
tons from the CERN-SPS impinging on a beryllium rod 
of 40 cm length and 2 mm diameter positioned 16 m up- 
stream of the experimental target. A magnetic channel 
selected negative particles with a mean momentum of 345 
GeV/c and a momentum spread of a(p)/p = 9%. A beam 
of 4.0 x 10 10 protons per 2.1 seconds spill in a machine 
cycle of 14.5 seconds yielded about 4.5 x 10 5 ir~ and 1.8 x 
10 5 Yi~ at the experimental target. A transition radiation 
detector was used to suppress ir~ on the trigger level j|. 
A detailed description of the hyperon beam setup and pa- 
rameters can be found in Q. 

The experimental target consisted of one 4 mm thick 
copper plate and three 2.2 mm thick carbon plates spaced 
at 2 cm intervals along the beam. For the carbon targets 
industrial diamond with a density of 3.3 g/cm 3 was used 
instead of graphite since it allows for thinner targets at the 
same interaction length, thus reducing the decay proba- 
bility inside the target. The total interaction length of all 
targets was 4.95%. The spacing of 2 cm ensured that most 
charm decays occurred in the gaps between the targets 
and hence had a lower background from secondary inter- 
actions. 

A total of 23 silicon microstrip detectors with a pitch 
of 25 and 50 /im, assembled closely downstream of the tar- 
gets, was used for the reconstruction of the charm produc- 
tion and decay vertices. Another set of 14 planes upstream 
of the targets was used for beam track reconstruction || . 

The vertex area was followed by a 12 m long decay zone 
for short-living strange particles, equipped with 38 drift 
chamber planes with a wire spacing of 5 cm, a sensitive 
area of 80 x 80 cm 2 and a spatial resolution of 300 /im. 
The mean efficiency was 88%. To increase the tracking 
efficiency in the central region with its high track density 
the drift chambers were interleaved with 5 sets of MWPCs 



3 



WA89 

Hyperon beam 
1993 layout 




DC,MWPC MWPC DC Hod 1 Hod 2 

Target area Decay area Omega spectrometer RICH Calorimeters 



Beam Si 



TRD 



50micron 

Scin 



Beam Si Targets 



Scin 



50micron 



Vertex Si 



Scin 

12 x 25micron 6. 5x6. 5cm 
12 x 25micron 5.0x5.0cm 
5 x 50micron 5.0x5.0cm 
Angles: 0,90,-45,+45 



Target area 

Fig. 1. Setup of experiment WA89 in the 1993 run. The lower part shows an expanded view of the target area. 



of 4 planes each. These MWPCs had a wire spacing of 
1 mm, a sensitive region of s» 12 x 12 cm 2 and an efficiency 
of 99%. 

To improve the track connection between the target 
and decay regions three sets of 4 MWPCs with a wire 
spacing of 1 mm were installed 2 m downstream of the 
target. 

Particle momenta were measured in the ^-spectrometer, 
a super-conducting magnet with an integrated field of 
7.5 Tm. A total of 45 MPWCs inside the magnet and 8 
driftchamber planes and 4 MWPCs at the magnet exit al- 
lowed a momentum resolution of a(p)/p 2 ~ 10 _4 (GeV/c) _1 . 

Charged particle identification was accomplished by a 
ring imaging Cherenkov counter with a 5 m long radiator 
volume filled with nitrogen at atmospheric pressure. The 
Cherenkov threshold was at 7 = 41. UV photons were de- 
tected in driftchambers filled with TMAE-saturated ethy- 
lene. A resolution of a = 2.8 mm and a mean number of 
15.5 photoelectrons per ring, equivalent to a quality factor 
of Nq = 53 cm^ 1 , allowed pion/kaon and pion/proton sep- 
aration up to a momentum of 90GeV/c and 150GeV/c, 
resp. with a rejection factor of 10 ||. 

An electro-magnetic lead glass calorimeter and a hadro- 
nic lead/scintillator "spaghetti" calorimeter were located 
downstream of the RICH. Both were not used in the anal- 
ysis presented here. 



The trigger was relatively open and had an average 
efficiency of 60-70% for the charm decays discussed here. 
It required a minimum track multiplicity in the target re- 
gion, derived from two scintillators placed after the first 
12 microstrip planes. Also two high-momentum particles 
were required with their momenta estimated from hit cor- 
relations in scintillators and wire chambers inside and be- 
hind the magnet. 

The experiment had two main data taking periods in 
1993 and 1994 during which a total of 350 million inter- 
action triggers were recorded, corresponding to an inte- 
grated luminosity of 450 ^b -1 on copper and 1160 /ib -1 
on carbon. 



3 Data Analysis 

Since there were only minor differences in the detector 
setups of 1993 and 1994, both data sets could be han- 
dled identically. The reconstruction of charm particle de- 
cays followed a candidate-driven approach: for each track 
combination in the target area whose charges matched 
the decay considered, reconstruction of the decay vertex 
was attempted. If the reconstruction was successful, the 
production vertex was reconstructed from the remaining 
tracks rejecting those which contributed most to the \ 2 
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Decay mode 


branching ratio 


events/background 


n 


b [c'/GeV 2 ] 


D u 


-> K-7T+ 


(3.83 ±0.12)% 


140/54 


6.2 ±0.8 






— > K~7r + 7r + 7r _ 


(7.5 ± 0.4)% 


86/29 


5.5 ±0.8 




DO 


-> K+7T 


(3.83 ±0.12)% 


195/40 


4.5 ±0.4 






— + K + 7T _ 7r + 7T _ 


(7.5 ± 0.4)% 


99/21 


4.6 ±0.6 




D+ 


— > K _ 7r + 7T + 


(9.1 ±0.6)% 


116/30 


5.4 ±0.6 


1.15 ±0.15 


D~ 


— > K + 7T _ 7T _ 


(9.1 ±0.6)% 


187/29 


4.2 ±0.4 


1.15 ±0.15 


D" 


K+K-tt" 


(4.6 ± 1.2)% 


55/15 


4.0 ±0.8 


0.79 ±0.12 


A c 


-> K-TT+p 


(4.4 ± 0.6)% 


49/15 


3.0 ± 1.0 


0.78 ± 0.20 



Table 1. Numbers of observed events and values of the production parameters n and b in the decay channels analyzed. Branching 
ratios are taken from M. 



of the vertex fit until x 2 / (degrees of freedom) fell below a 
given limit. 

Opposed to the total analysed sample which contains 
about 350 million events, the final charm samples show a 
total of only a few 100 entries per histogram. The reduc- 
tion factor of more than 10 6 was accomplished in three 
passes by cuts on (values in brackets give typical cut val- 
ues and resolutions): 

— the separation between production and decay vertex 
(> 7-13cr; crw 500 fxm), 

— the impact parameter of the reconstructed charm can- 
didate track w.r.t. the production vertex (< 3—6 a; a w 
lOjum), 

— the impact parameters of the decay vertex tracks w.r.t. 
the production vertex (> 3 — 6 er; a » 10 /im) 

Additional requirements were: 

— a decay vertex location outside the targets 

— a soft RICH identification of protons and kaons 

Table |l| lists the reconstructed decay modes, the num- 
ber of reconstructed events and the respective branching 
ratios. Figure || shows the corresponding signals. 

Acceptances and reconstruction efficiencies were de- 
termined using a detailed Monte Carlo simulation of the 
experiment. The charm events were simulated by an in- 
clusive event generator which generated the charm (anti- 
charm) particle considered according to the measured kine- 
matic distributions. The associated anti-charm (charm) 
particle was produced with the same kinematic distri- 
butions as the first particle, taking into account the az- 
imuthal correlation and 7"] p 2 distribution measured in a 
pion beam experiment @]7 For the associated particles a 
mixture of 25% of charged D-mesons, 50 % of neutral 
D-mesons and 25 % of A c was taken, in accordance with 
average yields of charmed particles measured in various 
experiments. The efficiencies were not sensitive to these 
percentages. The decay channel of the first particle was 
chosen explicitly while the decay of the associated charm 
particle was generated according to the known branching 
ratios. The remaining multiplicity of the event was gen- 
erated by simulating the interaction of one hadron carr 
ing the remaining momentum, using the Fritiof 7.02 
event generator for hadron-nucleus interactions. The de- 
tector response was generated using Geant 3.21. About 
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5 



T T — ' ' ▼ ■ — 1— 



E .r 



10 



■2-t- 



T acceptance 
* acceptance + trigger 
■ acceptance + cuts 



10 



0.2 0.4 0.6 0.8 1 



WlO 



10 



10 




▼ acceptance 

A acceptance + trigger 

■ acceptance + cuts 

• all 



pf [GeVV 2 ] 



Fig. 3. D detection efficiency as a function of xf and p t . "cuts" denotes the charm identification cuts. 
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Fig. 4. D differential production cross sections da/dxF and da /dpi (arbitrary units). Fits as described in the text. 



80000 events were simulated per channel with different xf 
distributions to ensure sufficient statistics at high xf- 

The efficiencies obtained for the different decay chan- 
nels analyzed show a similar behaviour. As an example, 
figure U shows the efficiencies for D~ as functions of pi and 
xf- While the efficiencies are practically flat in pf over 
the whole relevant region, the dependence on shows a 
sharp rise at about xf — O.f and levels off above xf = 0.2 
at a value of f 0% - 20%, depending on the decay channel. 
At high xf, the charm identification cuts produce a de- 
crease in efficiency because of the difficulty to reconstruct 
a production vertex from the few remaining low momen- 
tum tracks, which often are not seen in the detector. 

After the TRD decision in the trigger only « 2/3 of 
the beam particles are A large amount of contami- 
nation originates from pions from E~ decays between the 
end of the beam collimator and the target. These were 
suppressed by a cut on the track position/angle correla- 
tion which is fulfilled for beam particles coming from the 
production target, but not for decay pions (see Q for de- 



tails). Table g shows the beam composition after this cut 
for both data taking periods. The remaining contributions 
are fast pions, misidentified as X - in the TRD, and a few 
K~ and S~ which cannot be distinguished from S~. 



beam component 


1993 


1994 




0.14 ±0.012 


0.21 ±0.012 


K" 


0.038 ± 0.006 


0.048 ± 0.008 




0.013 ±0.001 


0.012 ±0.001 


E~ 


0.81 ±0.02 


0.73 ± 0.02 



Table 2. Beam composition after elimination of E decays. 



To correct for the beam contamination, we have to 
make assumptions about the dependence of the charm 
production cross section on the beam projectile, fn ex- 
periment E769 |l^] the ratio of D/D production by ir~ to 
production by protons was measured to be r vp = a(nN — * 
D/~D)/a(pN -> D/D) = 1.31 ± 0.24 at 250 GeV/c . With 
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a 360 GeV/c ir~ beam and a 400 GeV/c proton beam, the 
LEBC-EHS (NA27) collaboration measured this ratio to 
be r np — 1.05 ± 0.20 fyl. For the total charm production, 
a ratio o(ttN — > cc)/cr(pN — ► cc) = 1.00 ± 0.15 was ob- 
tained in NLO and NNLO QCD calculations We 
therefore correct for the pion and kaon contamination by 
using a value of rv P = 1.2±0.3, which amounts to a shift of 
(4.5±7)%. The possible effect from the S~ contamination 
is negligible in this context. 

4 Cross Section Results 

Since the detector acceptance drops sharply at small val- 
ues of xp, the measured differential cross sections were 
extrapolated to xf = using the usual ansatz da/dxF oc 
(1 — xf) 71 where n was obtained from a fit to the data 
above xf = 0.1. As an example, the measured da/dxp 
distribution for D~ production is shown in figure to- 
gether with the fit results. The values of n obtained in 
the fits are included in table |l[ The errors quoted are the 
statistical errors of the fits only. The values of n obtained 
are in agreement with values obtained in pion and proton 
beam experiments ( |T(],|l6| and earlier references therein). 
The errors preclude any conclusions about possible differ- 
ences in n between the different charm particles produced 
in this experiment or between our E~ beam result and the 
pion or proton beam results. 

We also fitted the observed pt spectra with a form 
do /dpi cx exp(—b ■ pi), assuming b to be independent of 
n. Figure |] shows the spectrum and fit result for D~ pro- 
duction, and the values of b obtained are listed in the last 
column of table [l]. Again we observe agreement with the 
earlier experiments quoted above, but are unable to draw 
further conclusions. 

The data from both targets were added and extrap- 
olated from S~ -nucleus to S" -nucleon cross sections 
assuming an A a dependence with a=l. The same depen- 
dence was used in the analyses of the proton beam exper- 
iments listed in table The most precise measurements 
of the exponent a in charm production come from the 
pion beam experiments WA82 @, E769 |J] and BEAT- 
RICE(WA92) The averages from the three experi- 
ments are a = 0.95 ± 0.04 for both D°/D° and D ± pro- 
duction, with no significant dependence of a on the lon- 
gitudinal or transverse momenta. Use of a = 0.95 would 
increase the resulting £~ -nucleon cross sections by a fac- 
tor of 1.13 (1.23) for carbon (copper) targets. 

The resulting integrated cross sections are listed in ta- 
ble ||. A first consistency check of the analysis can be ob- 
tained by comparing the results from the two-body and 
four-body decay modes of the D°, which are in good agree- 
ment. 

The total cc production cross section is taken to be 
the mean of the sums of the observed charm and an- 
ticharm cross sections (using a weighted mean for the 
two D° cross sections) with the result: <J C c(xp > 0) = 
5.6 ± 0.4 (stat) fib. 

Extensive studies of possible systematic effects have 
been carried out. No considerable dependence of the result 



Decay 


a(c)\ub] 


a(c)\ub] 


D u 


-> Ktt 


2.93 ± 0.29 


2.35 ± 0.33 




— > K7r7T7r 


2.07 ± 0.41 


2.38 ± 0.52 




mean 


2.64 ± 0.24 


2.36 ± 0.28 


D± 


— > K7T7T 


1.30 ± 0.14 


1.13 ± 0.16 


D± 


-> KKn 


1.60 ± 0.42 




A c 


— > Kp7T 




2.37 ± 0.68 


Sum 


5.54 ± 0.51 


5.86 ± 0.75 



Table 3. Production cross sections for anti-charmed and 
charmed hadrons in the range if > . The third line gives 
the D° cross sections averaged from the measurements in the 
K-7T and Y^-ktt-k decay channels. Errors are statistical only. 



on the cut values or on reasonable variations of detector 
efficiencies was found. We assign an overall systematic er- 
ror of 20% to all cross sections due to the uncertainties 
of the Monte Carlo simulation, the extrapolation to xf=0 
and the luminosity measurement. 

Our cross section result is based on the observation 
of D and D s mesons, to which all excited charm mesons 
will eventually decay, and of A c , to which all charmed 
non-strange baryons will decay. Therefore possible unde- 
termined contributions from charm anti-baryons and from 
charmed-strange baryons remain, which will be discussed 
in the following. 

A significant contribution to the total cross section by 
charm anti-baryons is not expected due to the observed 
strong asymmetry in A c /A c production ||. From the re- 
sults of a direct search for charmed-strange baryon pro- 
duction, we cannot exclude production cross sections of 
1 fib |17]] . We therefore add a further systematic uncer- 
tainty of +1 fib from S c production to our total cc pro- 
duction cross section. A contribution from f2 c should be 
small compared to the overall uncertainties. 
The total cc production cross section for xf > then is: 
a c7 :(x F > 0) = 5.6 ± 0.4 (stat) ± 1.1 (syst) + 1.0 (H c ) fib 

5 Discussion 

In the hyperon beam experiment WA89 at CERN, close to 
1000 charm particles were found in 340 GeV/c S _ -nucleus 
interactions. We determined production cross sections for 
D°, 130, D*, D- and A+ in the region x F > 0. While a 
clean A+ signal could be observed, the A c could not be 
singled out. This suppression of antibaryon production in 
baryon beams is well known from non-charm baryon pro- 
duction. The cross sections for and D/D production 
are comparable. The suppression of Dj~ production rel- 
ative to DJ is probably caused by the suppression of s 
production relative to u, d production, while the s quark 
of the DJ can be supplied by the E~ beam particle. From 
this argument one would expect an enhancement of D~ 
over D + production as well, but this enhancement is not 
so strong since light antiquarks are produced more co- 
piously than the heavier s quarks. Indeed, we observed a 
significant D _ /D + production asymmetry at larger xf [§, 
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Experiment 


beam, 
mom. 


targets 


a(D/D) 
[/xb] 


xp range 


<t(cc) for xp > 
[Mb] 


NA32 


IN 


p, 200 


Si 


1.5±0.7 


>0 


2.25±1.10 


E769 [ 


LO 


h 


p, 250 


Be,Al,Cu,W 


9.0±1.5 


>0 


13.5±2.9 


NA27 


[L 




p, 400 


P 


30.2±2.9 


all 


22.6±3.8 


E653 [ 


L9 




p, 800 


emulsion 


76±20 


all 


57±17 


E743 [ 


20 




p, 800 


P 


48±12 


all 


36±10 


WA89 


E" , 340 


C,Cu 


7.4±1.5 


>0 


11.3±2.4+2.0 



Table 4. Production cross sections measured in proton beams and in this expt. Beam momenta are given in GeV/c . 
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Fig. 5. D/D production cross sections by E~ (this expt., full 
circle) and protons (open circles). 



but the effect on the integrated cross section for xf > 
is not observable within the errors of our measurements. 

In figure || and in table f| our result for the (D or D) 
production cross section a B ^(xF > 0) = 7.4±0.4(stai.)± 
l.bj syst.) [ib is compared to measurements in proton beams 

EM 
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Proton beam measurements valid for — 1 < 
Xp < 1 have been divided by 2 for this comparison. Our 
result is significantly lower than the proton beam values. 
It is not clear, however, whether for this comparison we 
should include all or part of our measured value for D~ 
production, a D - = 1.6 ± 0A(stat.) ± 0.3(syst.) /j,b, since 
our beam projectile carries a strange quark. 

Our measured value for the ratio of charged and neu- 
tral D/D production, cr(D ± )/cr(D°, D°) = 0.49 ± 0.06, is 
compatible with the mean value from the proton beam ex- 
periments gapMpO), cr(D ± )/cr(D°,D0) = 0.66 ±0.09. 



The corresponding ratio from pion beam experiments at 
beam energies between 200 and 600 GeV/c is cr(D±)/ 
er(D°,D°) = 0.415 ± 0.010 in good agreement with our 
result (for the pion beam results, see ref. |n6[ and ear- 



lier references therein). For the ratio of D s and D pro- 
duction we obtained a(D~)/a(D,T5) = 0.22 ± 0.06. This 
ratio has not been measured in proton beams. In the 
pion beam experiments quoted above, the ratio of D+ and 
D~ production to D production has been measured to be 
er(D±)/cr(D,D) = 0.129 ± 0.012. Since both D+ and D" 
are non-leading when produced in a pion beam, we can 
safely assume that <t(D~)/ct(D, D) = 0.065 ±0.006 in pion 
beams, which is markedly lower than our S _ beam result. 
This is another clear manifestation of the leading particle 
effect. 

Theoretical calculations exist only for the total charm 
production cross section, cr C c, for all x^. Therefore our re- 
sult and the results of the proton beam experiments | fio| 
pl| , obtained in the kinematic region xf > 0, were mul- 
tiplied by a factor of 2. This factor of 2 holds for inter- 
actions with equal beam and target particles, but even 
in X - N interactions, no significant deviation from 2 is 
expected, since the production cross sections integrated 
over the range < Xp < 1 show only small differences 
between leading and non-leading charm hadrons as dis- 
cussed above. Furthermore, the D/D production cross sec- 
tions measured in the proton beam experiments have to 
be extrapolated to total charm production. Our value for 
the ratio of total charm production to D/D production is 
cr c/ c/cr D/S = 1.53±0.15+0.10. The E769 collaboration go) 

measured o" d+ j5 = 9.0 ± 1.5^(6 and 5.1 < <j\ c < 21.5/i6 in 
their proton beam, which is consistent with our ratio. In 
their pion beam they measured er(A c + D s ) / cr(D + + D°) = 
0.53±0.14 (particles and antiparticlcs added). In the pho- 
toproduction experiment NA14/2 this ratio was measured 
to be 0.36±0.11 pl| . These values indicate that the ratio 
cr c / 5 /<7 D ,jj does not depend strongly on the nature of the 

beam particle. We used a value o' c /c/ (J d/'d = ^ 0-2 1° 
calculate a c /c from the proton beam results on <r D /g. Fi- 
nally, the production cross section for cc-pairs, a C c, was 
obtained by dividing the above results for c or c by a fac- 
tor of 2. The results are listed in the last column of table ^ 
and are compared in figure |^ with theoretical calculations. 
The dashed/solid/dotted bands were obtained by Frixione 
et al. in NLO QCD-calculations for charm quark mass val- 
ues of 1.2/1.5/1.8 GeV/c 2 , respectively |§. The width of 
the bands for fixed m c is given mainly by the uncertainties 
of the renormalization scales. The proton beam measure- 
ments prefer the upper range of the predictions, and seem 
to exclude a value of m c — 1.8 GeV/c 2 . This holds also 
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Beam momentum (GeV/c) 



Fig. 6. Total charm production cross sections by E~ (this 
expt., full circle) and protons (open circles) compared to the- 
oretical predictions as explained in the text. 



for charm pro duction measurements in pion beams, as dis- 
cussed in |12|, 22 1 . The thick dot-dashed line denotes the 
results of Smith and Vogt [|l3| for a charm quark mass of 
1.5 GeV/c 2 . These authors used resummation techniques 
and approximations for NNLO. At our beam energy, their 
results are larger by a factor 3 than the NLO calculations 
for the same charm quark mass. 

Our experiment provides the first comparison of charm 
production by £~ and by protons. Our result seems to 
be lower than the proton beam results, in line with the 
observation that the total S~ N cross sections are smaller 
than the pp cross sections [E3|. However, in view of the 
problems encountered when comparing cross sections from 
different experiments, it may be premature to claim that 
er C c or & D /~d is smaller in £~ than in proton beams. The 
forthcoming results from the hyperon beam experiment 
SELEX at FNAL will help to clarify this point & 



References 

1. P. Nason, S. Dawson, R.K. Ellis, Nucl. Phys. B 327, 49 
(1989) 

2. WA89 collaboration, M.I. Adamovich et al., Eur. Phys. J. 
C 8, 593 (1999) 

3. W. Bruckner et al., Nucl. Instr. Meth. A 378, 451 (1996) 

4. Yu.A. Alexandrov et al., Nucl. Inst. Meth. A 408, 359 
(1998) 

5. W. Bruckner et al., Nucl. Instr. Meth. A 348, 444 (1994), 
Nucl. Instr. Meth. A 357, 274 (1995) 

6. W. Beusch et al., Nucl. Instr. Meth. A 323, 373 (1992); 
H.-W. Siebert et al., Nucl. Instr. Meth. A 343, 60 (1994); 
U. Miiller et al., Nucl. Instr. Meth. A 371, 27 (1996) 

7. The particle data group, Phys. Rev. D 54, 1 (1996) 

8. BEATRICE collaboration, M.I. Adamovich et al., Phys. 
Lett. B 385, 487 (1996) 

H. Pi, Comput. Phys. Comm. 71, 173 (1992) 
E769 collaboration, G.A. Alves et al., Phys. Rev. Lett. 77, 
2388 and 2392 (1996) 

LEBC-EHS collaboration, M. Aguilar-Benitez et al., Z. 
Phys. C 40, 321 (1988) 
S. Frixione et al., CERN-TH/97-16 
J. Smith and R. Vogt, Z. Phys. C 757~271 
WA82 collaboration, M.I. Adamovich et al., Phys. Lett. B 
284, 453 (1992) 

E769 collaboration, G.A. Alves et al., Phys. Rev. Lett. 70, 
722 (1993) 

BEATRICE collaboration, M.I. Adamovich et al., Nucl. 
Phys. B 495, 3 (1997) 



9. 
10. 

11. 

12. 
13. 
14. 

15. 

16. 

17. 

18. 

19. 

20. 

21. 
22. 
23. 
24. 



hep-ph/9702287 
(1997) 



C. Busch, thesis, Heidelberg 1998 (unpublished); H. Neeb, 
thesis, Heidelberg 1998 (unpublished) 
ACCMOR collaboration, S. Barlag et al., Z. Phys. C 39, 
451 (1988) 

E653 collaboration, K. Kodama et al., Phys. Lett. B 263, 
573 (1991) 

LEBC-MPS collaboration, R. Ammar et al., Phys. Rev. 
Lett. 61, 2185 (1988) 

M.P. Alvarez et al., Z. Phys. C 60, 53 (1991) 

S. Frixione et al., Nucl. Phys. B 431, 453 (1994) 

S.F. Biagi et al., Nucl.Phys. B 186, 1 (1981) 

SELEX collaboration, J. Russ et al., CMU-HEP98-07 and 



hep-ex 9812031 



